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Abstract 

The key quantity of the heavy quark theory is the quark mass mg. Since quarks are 
unobservable one can suggest different definitions of mg. One of the most popular choices 
is the pole quark mass routinely used in perturbative calculations and in some analyses 
based on heavy quark expansions. We show that no precise definition of the pole mass can 
be given in the full theory once non-perturbative effects are included. Any definition of this 
quantity suffers from an intrinsic uncertainty of order AgcD/ m Q- This fact is succinctly 
described by the existence of an infrared renormalon generating a factorial divergence in 
the high-order coefficients of the a s series; the corresponding singularity in the Borel plane 
is situated at 2n/b. A peculiar feature is that this renormalon is not associated with the 
matrix element of a local operator. The difference A = Mh q — mg° le can still be defined in 
Heavy Quark Effective Theory, but only at the price of introducing an explicit dependence 
on a normalization point fi: A(/x). Fortunately the pole mass toq(O) per se does not appear 
in calculable observable quantities. 
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1 Introduction 



Significant progress has been achieved recently in the quantitative treatment of the decays 
of heavy flavor hadrons by employing expansions in powers of 1/mn, where mq denotes 
the heavy quark mass. Exclusive transitions between two heavy flavor hadrons are conve- 
niently dealt with by using the spin- flavor symmetry of heavy quarks Q (see also ref. [0]); 
the formalism of Heavy Quark Effective Theory (HQET) ||, ||] incorporates this symme- 
try concisely at the Lagrangian level. Inclusive decays, on the other hand - nonleptonic, 
semileptonic or radiative transitions to any type of the final state - can be treated directly 
in QCD [||, ||, fj], ||, |9], [K| via Wilson's operator product expansion 11]. 

The key parameter in most aspects of heavy quark physics is obviously 'the heavy 
quark mass'. There are no problems in defining this quantity within purely perturbative 
calculations. At this level one choice turns out to be particularly convenient: it is the 
so-called pole mass of the heavy quark, defined as the position of the pole in the quark 
propagator in perturbation theory. This quantity, introduced in QCD in the seventies (see, 
e.g. Ref. [p^| ), is well defined in each finite order of perturbation theory; unlike many other 
definitions, it is introduced in a gauge invariant way. This convenient feature has made it 
very popular. Important results of perturbative calculations such as the total semileptonic 
widths (including radiative corrections) are routinely expressed in terms of the pole mass 
(see e.g. Q). 

In this paper we exhibit an important drawback in the concept of a pole mass that 
becomes apparent as soon as one addresses leading non-perturbative corrections to order 
1/rriQ. The problem arises, of course, because the pole mass is sensitive to large distance 
dynamics - although this fact is not obviously seen in the standard perturbative calculations 
- and the corresponding treatment requires special care. It had actually been noted before 
that marrying full QCD with the notion of the pole mass faces subtle difficulties (TJ . 

Our main assertion will be three-fold: 

(A) Infrared contributions lead to an intrinsic uncertainty in the pole mass of order 
Aqcd, i-e. an effect of relative weight 1/mg. Perturbation theory itself produces clear 
evidence for this non-perturbative correction to m P Q le . The signal is the peculiar factorial 
growth of the high order terms in the a s expansion corresponding to a renormalon [15, |l^] 
residing at 2ir/b in the Borel plane, where b is the first coefficient in the Gell-Mann-Low 
function, b = (UN c /3) — (2Nf /3). The physical reason lying behind this linear effect is just 
the "Coulomb" energy of the heavy quark. 

A subtle point is that infrared-renormalon effects can usually be associated with the expec- 
tation values of some local operators [JO], [lTj of the corresponding dimension. This general 
pattern is not realized in the case of the pole mass for the reasons which will be clarified 
below. 

(B) The Heavy Quark Expansion (HQE) yields for directly observable quantities such 
as the total semi-leptonic width of the heavy flavor hadron decay: 
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where for definiteness we consider B meson decay and the coefficient c is c = 1 + 0(a s ). In 
the total semi-leptonic width non-perturbative effects of order 1/toq are absent (see Refs. 
[0, [|), and the corrections start from l/rng (the factor (2Mb)" 1 reflects the relativistic 
normalization of the state | -£>}). While the pole mass is useful within purely perturbative 
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calculations, it makes no sense to employ it in such an analysis that aims for an accuracy 
~ 0(Aqcd/itT'q) or even ~ ©(Ag^/mg) - simply because it cannot be unambiguously 
defined at order Aqcd / m Q- 

(C) There is a profound theoretical difference between the pole mass and the total 
semileptonic width: calculation of the latter can be formulated as an OPE-based procedure, 
ensuring factorization of the large and small distance contributions; the pole mass, however, 
cannot be treated in such a way. 

On the other hand - and fortunately - there is no need to use the pole mass in these 
inclusive transition rates. Careful analysis shows that contrary to popular opinion it is the 
running mass m(fi) with /j, 3> Aqcd that naturally enters. The properly defined running 
mass includes only the effects of momenta higher than /i; therefore if /i is chosen sufficiently 
high there is no infrared uncertainty in m(^). In particular, the parameter most relevant 
to the inclusive heavy quark decays is m{mo). In other problems it may turn out that the 
running mass normalized at a different point enters; each particular case requires its own 
careful analysis. It is important however, that this normalization point never goes down 
to a typical hadronic scale if one wishes to avoid the corresponding infrared uncertainty in 
the Wilson coefficients. We stress here that the normalization point fi does not necessarily 
coincide with the off-shellness of the heavy quark inside the heavy hadron. 
Certainly, the definition of the running mass m(/i) depends both on the scheme and the 
gauge used - an obvious inconvenience. In this respect, however, the running mass does not 
differ - as a matter of principle - from the running gauge coupling constant a s (/i) where an 
explicit scheme and gauge dependence first emerges on the third loop level. The essential 
difference is that for m(/x) this dependence manifests itself already at the one-loop level; 
the observables, however, are independent. 

Correspondingly, any consistent definition of the parameter A must explicitly introduce the 
normalization point /i, so that A is actually running, A(//), although the running law is 
somewhat unconventional, see eqs. ji~7|), (]48|). 

The remainder of this paper will be organized as follows: in Sect. 2 we restate the general 
procedure of separating off the infrared effects within Wilson's OPE and the relation between 
the infrared renormalons and non-perturbative condensates; in Sect. 3 the perturbative 
series for the heavy quark mass is discussed; Sect. 4 is devoted to the infrared renormalon 
contribution to the mass; Sect. 5 demonstrates that the pole mass does not appear in 
directly observable quantities, such as, say, the total widths; finally in Sect. 6 we discuss 
the running of A and summarize our results. 



2 Separation of the Infrared Effects within Wilson's OPE 

In this section we remind the reader of a crucial property of the operator product expansion 
[ |l7| , 18, IS]: it provides us with a systematic separation of infrared and ultraviolet contri- 
butions. This brief excursion into standard OPE applications will allow us to reveal a basic 
difference in how infrared effects enter in the pole mass and, say, into correlation functions 
at large momentum transfer. 

As an example let us consider the correlation function of vector currents, 

V = {q^u ~ «V) n = i J d 4 xe^(0|T{j M (x)>(0)}|0> (2) 
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where 

with ip being a massless quark field. To avoid "external" logarithms which are irrelevant 
for the problem under discussion one usually deals with a modified quantity defined as 

fl = -Air 2 Q 2 (d/dQ 2 )U , Q 2 = -q 2 . (3) 

The one- and two- loop graphs determining II in perturbation theory are depicted in Fig. 1. 
The well-known calculation of these diagrams yields 

ft-l + 2!<92. (4) 

The virtual momenta saturating the corresponding loop integrals are of order Q, and if Q 2 
is chosen to be large the characteristic virtual momenta are also large. The result (||) then 
represents the short distance contribution. 

The fact that eq. @) is correct in perturbative QCD does not mean, however, that it is 
correct in the full theory. Indeed, in deriving eq. (||) one integrates over all gluon momenta 
k, including the domain of small k where k is the momentum flowing through the gluon 
line in Fig. 16. Of course, this domain gives a relatively small contribution to the integral; 
nevertheless this contribution is definitely wrong since for small k the Green functions are 
strongly modified by non-perturbative effects and have nothing to do with the perturbative 
propagators one uses in obtaining eq. (ffl). 

The emergence of non-perturbative corrections can actually be inferred from pertur- 
bation theory per se if one recalls the presence of the Landau singularity in the running 
coupling constant in the infrared domain (for QCD). This is the essence of the concept of 
infrared renormalons [15, 16]. The phenomenon is quite simple and manifests itself in the 



behavior of the high order terms in the a s expansion. The relevant diagrams are those 
where a chain of loops has been inserted into the gluon Green function (Fig. 2). This chain 
is equivalent to replacing the fixed coupling by the running coupling constant a s (k 2 ) in 
the integrand and integrating it over with some weight function; the latter is given by the 
remaining propagators in the diagram. This weight function is such that the integral over 
d k converges in the ultraviolet domain. To illustrate how this works let us use, following 



Ref. [20 1, a simplified expression for II, 



fr n2 [ju2 k 2 a s (k 2 ) 

LL renorm ~ V J UK ^ g . [0) 

This expression coincides with the original one in the limits k 2 <C Q 2 and k 2 ^> Q 2 relevant 
to the infrared and ultraviolet renormalons, respectively. 

The integral in eq. @ is saturated at k 2 ~ Q 2 ; to get the main contribution one 
substitutes a s (Q 2 ) for a s (k 2 ) and finds the standard two-loop result for ft. If, however, we 
are interested in high orders in a s (Q 2 ) the tails at k 2 <C Q 2 and k 2 ^> Q 2 become important. 
The contribution of these tails have the form 

rQ 2 

U IR = Q- 4 k 2 dk 2 a s (k 2 ), (6) 

JO 

and 
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where the subscripts are self-evident. Substituting the running gauge coupling 



a s (k 2 



l-(&a s (Q 2 )/4^)ln(Q 2 /fc 2 )' 



where b is the first coefficient in the Gell-Mann-Low function, and expanding in a s (Q 2 ) we 
get the whole series in a s (Q 2 ). Namely, 




(8) 



and 
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(9) 

The first expression corresponds to the infrared renormalon while the second one refers to 
the ultraviolet renormalon. The factorial growth of the coefficients is explicit in both cases. 
There are two differences, however. The ultraviolet renormalon is represented by the sign- 
alternating series, in distinction to Hir. The second difference is the n dependence of the 
coefficients in front of n!(6ct s /47r) n+1 . In the infrared renormalon this coefficient is 2~^ n+l \ 
compared to 1 in the ultraviolet renormalon. These differences result in the fact that the 
positions of the singularities in the Borel plane (to be denoted b) are at 

~ 8-7T ~ 47T . . 

In what follows we will not consider the ultraviolet renormalon, as well as other sources 
of the singularities in the Borel plane, for instance, instanton-antiinstanton pairs (see e.g. 
producing a singularity at b = Air. 

The occurrence of the factorial in Tim is correlated with the fact that at large n the 
integral is saturated not at the large scale k 2 ~ Q 2 but, instead, at parametrically low scales 
k 2 ~ Q 2 /n. Therefore when n becomes large one encounters the Landau singularity, and - 
in an indirect way - comes to the realization that in the infrared domain the gluon Green 
function cannot coincide with the perturbative expression. 

The factorial growth of the coefficients points to the fact that one cannot - as a matter of 
principle - infinitely increase the accuracy of the perturbative approximation by including 
higher and higher terms in the series: starting from some order the absolute size of the 
corrections increases again. The best one can achieve is to truncate the series at an optimal 
value of n ensuring the best possible approximation. It is not difficult to check that the 
error introduced by this truncation is of the order of 

r 87T h-QCD 

Since the series is not Borel-summable, it is in principle impossible to achieve a better 
accuracy within the framework of perturbation theory. All one can conclude is that there 
must be a non-perturbative effect of the form (p4|). 

It is well known that this problem does not mean that one has to abandon hope for 
establishing theoretical control over the terms of this order of magnitude. The way out is 
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provided by a consistent use of Wilson's OPE procedure. First one introduces a momentum 
scale [i such that momenta above \i can be treated perturbatively. All loop momenta are 
classified according to whether they exceed \i or lie below it. The integration over the 
infrared domain below fi is then explicitly excluded from the perturbative calculation and 
one finds 

U pert = U\ k>fl = H Const -j . (12) 

7T 

The fact that the subtracted term in this particular case is proportional to /i 4 is not ac- 
cidental, of course: it can be anticipated from eq. (11) and it will be clarified shortly. 



This result automatically follows from the explicit calculation, provided that the scale [i is 
introduced without breaking the gauge invariance of the theory (in practice however this 
may turn out to be a technically highly non-trivial exercise). Excluding the domain k < /j, 
from the perturbative calculation does not mean that we just lose this contribution. Within 
the Wilson OPE the contribution of this infrared domain re-enters through the vacuum 
expectation value 

AU = -^(a s G 2 ) (13) 

where G = G® is the gluon field strength tensor, and the operator G 2 in the right-hand 
side is normalized at \i (i.e. by definition this vacuum expectation value includes all virtual 
momenta below The vacuum expectation value (a s G 2 ) has the form 

(a s G 2 ) = Const Aq CD + Const a s (/x 2 )/i 4 (14) 

and in the full expression H per t + All the dependence on the auxiliary parameter \x cancels. 
The requirement that the dependence on \i cancels in the end dictates that fi must appear 
in tip er t as fx , since no gauge invariant operator of lower dimension exists. 

Let us emphasize that the normalization point /j, should be high enough to ensure the 
applicability of perturbative calculations above /i, i.e. 

a 3 (fj) < 1. 

On the other hand, it is desirable to have this parameter as small as possible, so that 
the corresponding terms will represent insignificant corrections in eqs. (12), (|i~4T). If this 



wish can be satisfied there is no need to carefully work out details of how to introduce /j, 
explicitly! In particular, in eq. (p^|) the term with A.% CD should be much larger than that 
with /j, 4 . In other words, the numerical values of the condensates are assumed to be much 
larger than their perturbative parts for some [i belonging to the window discussed above. 
This is what is called the practical version of OPE - powers of fi do not show up explicitly 
then, although, of course, the conceptual necessity of having /i should be always kept in 
mind (/u is certainly kept in logarithms in the practical version of OPE). It is fortunate for 
applications of QCD that the practical version of OPE works well in the vast majority of 
instances; this could not have been anticipated a priori. Otherwise, all calculations based 
on OPE would be much less useful since it would be mandatory to explicitly construct the 
procedure of introducing \i. 

What then happens with the infrared renormalon within Wilson's procedure? If one 
defines the perturbative part with the infrared cut off at the point /i, the factorial growth 
of the coefficients in the perturbative series in a s stops starting from some value of n, 
n ~ Q 2 /fi 2 , since the integral "wants" to be saturated in the domain k 2 < fi 2 and this 



6 



domain is now simply eliminated from the perturbative sector. The price one has to pay is 
the introduction of a new, non-perturbative parameter, the gluon condensate. Once it is in- 
troduced, however, the perturbation theory is amended, the effects oc Q~ A become tractable 
and this accuracy is legitimate. If necessary, one goes a step further. Of course, corrections 
of higher order in Q~ 2 require the introduction of new higher-dimensional condensates. 

Let us emphasize that the actual uncertainty due to the infrared renormalon, eq. (JTT]) , 
should by no means be equated with the contribution from the gluon condensate. It is true 
that they are of the same order in the parameter Aqcd/Q and the infrared renormalon 
anticipates the appearance of the gluon condensate; the gluon condensate contribution, 
however, is much larger numerically, and this is the reason why the practical version of 
OPE is so successful. One can interpret this fact in terms of extremely strong distortions of 
the Green functions in the infrared domain [17]. In other words the modifications (compared 
to smoothly extrapolated perturbative Green functions) are not just of the order of unity 
but are much larger numerically. 

The discussion above makes it clear that in principle perturbation theory already signals 
the emergence of a non-perturbative correction in n of order Aq^/Q 4 . A hint is provided 
by the simplest diagram of Fig. lb. If we introduce \x and calculate the graph of Fig. lb 
according to the Wilson procedure we discover a correction a s /i 4 . In combination with the 
general fact of strong distortion of the Green functions in the infrared domain this a s /i 
correction triggers the rest of the machinery which eventually leads to the non-perturbative 
gluon condensate contribution in II. 

The lessons one can draw from this rather standard procedure [17, 18, [l{J are quite 
evident: although it is impossible to actually calculate non-perturbative contributions by 
analyzing the behavior of perturbation theory for a particular quantity, the fact that non- 
perturbative contributions exist and their particular form can be inferred. Moreover, one 
can find out what kind of non-perturbative terms are to be expected in a two-fold way: by 
studying low order perturbative graphs within the Wilson procedure and by inspecting high 
orders of perturbation theory as they are generated by the infrared renormalon. 

Below we will use the both lines of reasoning to show that the "pole mass" of heavy 
quarks - if treated in the context of problems where we intend to include non-perturbative 
effects - contains a piece of the order of Aqcd which must be considered as an intrinsic 
uncertainty. Unlike the standard case, however, this piece is not related to any matrix 
element of a local gauge invariant operator. Therefore, one cannot amend the perturbation 
theory based on m P Q le in the manner it is usually done. This happens due to the fact that 
the notion of the "pole mass" by itself is ill-defined. 



3 Perturbative Corrections to tuq 

Following the general strategy outlined above we start our analysis of the pole mass with 
the simplest perturbative graph of Fig. 3. Instead of a straightforward calculation of this 
graph [] we follow the Wilson procedure and introduce, first of all, the normalization point 
H such that the domain of virtual momenta k < fi is discarded. Thus, we are going to 
calculate an analog of fl per t in eq. (|T^) . 

1 The corresponding computation can, of course, be found in any textbook on quantum electrodynamics 
(QED); to this order there is no difference between QED and QCD, up to a trivial overall factor. 
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As already mentioned, separating the infrared and ultraviolet domains would in gen- 
eral require a rather sophisticated machinery Fortunately, the situation simplifies in this 
particular case because the diagram of Fig. 3 is the same as in QED. In QED we can just 
introduce a "photon mass" A, which, on the one hand, preserves Ward identities associated 
with current conservation, and on the other, suppresses the contribution of all virtual mo- 
menta below A. In this way, the infrared domain is automatically discarded. The photon 
mass A is to be identified with the normalization point /Li, cf. eq. (|l^). The procedure of 
introducing \x suggested here cannot be extended to higher loops. It is quite satisfactory, 
however, for our more limited purpose, namely to establish the presence of a correction of 
order I/thq in the pole mass. 

Accordingly we use the following expression for the gluon propagator 

<« = -^-f^)^^. (15) 

Here £ is the gauge fixing parameter and the term Mq/(Mq — k 2 ) ensures ultraviolet regu- 
larization (Mo is the ultraviolet cut-off). We have introduced it to be safe in the ultraviolet 
regime, although it will play no role in what follows. Since we are interested in the infrared 
domain we further assume that 

Explicit evaluation of the diagram of Fig. 3 yields 



= mg> | 1 + J0»^| + Const.) ) - |^A. (16) 
m Q 



Here niq is the bare mass (the mass at Mo), while the quantity on the left-hand side is the 
pole mass which depends on A. We remind the reader that the domain of virtual momenta 
k < A is absent in toq(A), and this leads to the dependence on A. 

Next, we identify A and /x, as explained above, and express the mass at one normalization 
point (we will refer to this mass as the running mass) in terms of that corresponding to 
some "starting" normalization point no, 

m Q {n) = m Q (iJ,o) + — (fio-n). (17) 

3 7T 



What is important in eq. (|17| ) is the occurrence of a correction of order I/tuq relative to 
the leading term. 

The procedure outlined above is not unambiguous in both its elements - the definition 
of the running mass and the specific manner in which the normalization point has been 
introduced. In principle, one can use other prescriptions. Let us mention, for instance, the 
suggestion of Ref. |2^] where the running mass m<g(/i) was introduced through a certain 
integral over the cross section for the process 7 + Q — > g + Q. The normalization point \x 
then enters as an upper limit in this integral (see |^] for more details; we plan to discuss 
this approach in a forthcoming publication). Within this procedure one gets an analog of 



eq. (17) which still contains a linear correction, albeit with a different numerical coefficient. 



2 As will be seen below, naive dimensional regularization cannot serve this purpose. 
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While the numerical value of the 1/wiq correction is thus scheme-dependent, its presence 
is not! 

The physical meaning of the 1/mn contribution is quite transparent. It is nothing else 



than the classical Coulomb self-energy of a static color source, see e.g. Ref. JT4 ] . The 
energy of the Coulomb field at a distance rg is given by 

2a s 1 

Ecoui = -=*■-. 18 
3 r 



Identifying ro with l//x we recover eq. (17). 



The Coulomb contribution to the mass can readily be derived directly from the graph 
in Fig. 3. In the limit of large rriQ (compared to the gluon virtual momentum k), i.e. in 
the static limit, the expression for the mass shift takes the form 

5m Q ^m Q -m Q ^) = -i- 9s J {k(j +ie){k 2 +ie) ( 19 ) 

with an ultraviolet cut-off in the k integration implied at /x. (The integral for SrriQ becomes 
linearly divergent in the static limit.) It is easily seen that the only surviving contribution 
in eq. ( |l9|) comes from l/(fco + ie) in the form of a term —i7r5(ko) . The fact that the 
static limit implies the vanishing of ko is quite evident by itself; the occurrence of the linear 
infrared effect under consideration can be traced back to this feature of the static limit. 

Performing first the integration over ko (which reduces merely to putting ko = 0) we 
arrive at 

8tt f d 3 k 1 
5m Q = -ja s J j-^ T2 . (20) 

If an ultraviolet cut-off is introduced through a factor /Xo/(Mo + ^- ) m the integrand we get 

5m Q = m P Q le - m Q (fi ) = ^a s // , (21) 

cf. eq. (|i~7|) where, to get the pole mass on the left-hand side, one should set [i = 0. Were 
the cut-off introduced in a "hard" way we would get a different coefficient in front of the 
linear term, of course, but the very fact of its presence would remain intact. 

The appearance of the term linear in ^ in eq. ( |l7| ) tells one, according to the discussion 
of Sect. 2, that there is a renormalon singularity in the perturbative series giving rise to a 
relative uncertainty in the pole mass of order 1/mg. Below we will demonstrate it explicitly. 

To conclude this section an important remark of a conceptual nature is in order. The 
calculation carried out above clearly reveals an important fact: while the OPE-like proce- 



dure routinely used in HQET resembles closely Wilson's OPE procedure [11], it has one 
very distinct feature. In the static limit all energy transfers to the heavy quark line vanish, 
implying that the time separations are always large. Physically that is quite transparent: 
the heavy quark after being placed at the origin as a static color source at t = — oo remains 
there at rest till t = +oo. The OPE procedure in the effective low-energy theory is then 
based on a separation of small spatial from large spatial momenta. Therefore, below tuq we 
actually deal with a three-dimensional version of OPE, which results in peculiarities that 
might seem strange, at first sight, to those who got used to the standard features of the 
four-dimensional OPE. 
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4 Infrared Renormalon 



Next one examines the impact of the high order corrections in a s to the pole mass. We 
again study the chain of loops inserted into the gluon propagator, this time in the graph of 
Fig. 3, see Fig. 4. Summing all these "bubbles" amounts to replacing a s in eq. (|20| ) by the 
running coupling a s (k 2 ) in the integrand, 



pole , n ott f d°k a s (k z ) 
The running gauge coupling is given by 

1 - (ba s //g /4tt lnGug/F ^ \ 4tt fc 2 / 



Substituting the expansion of eq. ( |23| ) into eq. ( p2[ ) we immediately obtain the following 
series: 

pole , s 4q s (/x ) ^ ^ ba s (/i ) ^ w 

m Q -™q(h>)= 3n Vol^C n ^ ^ j (24) 

where the coefficients C n are given by 

^-L 1 <25) 

At large n these coefficients grow factorially, 

C n = 2 n n\ . (26) 
In other words one can say that the position of the nearest singularity in the Borel plane 



16 1 is at b = 2ir/b. This series is not Borel-summable due to the presence of an infrared 



n ° ~ ( 27 ) 



renormalon; truncating it at the optimal value of no 

2tt 
ba s (nl) 

one arrives at an estimate of the irreducible uncertainty in m P Q le — mgf/Jo), 

A(m^ - m Q G* )) ~ ±n> « x P{-^y} ~ | A ^d- (28) 

Thus, we see that perturbative QCD indeed does not allow one to define m^ le with an 
accuracy better than Aqcd, i-e. an infrared effect linear in Aqcd / m Q- 

A simple way to explain the estimate pSj ) is to turn to the original integral (22) 

m Q l£ ~ m Q(W)) = ^a s (fj, ) f — ^- -j. (29) 

v Sir Jq i QQa(w) !_ Mn 

Previously we have just expanded the denominator in the powers of a s (uo), obtaining in 
this way the factorial behavior of the expansion coefficients related to the existence of the 
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pole singularity in the integrand. Now, instead, we regularize the singularity, say, by adding 
it in the denominator. Then eq. ( ^9|) acquires the imaginary part, 



Im (m P Q le - m Q (/i )) = ^ A QCD (30) 

where Aqcd parametrizes the position of the infrared pole in the running gauge coupling. 
The estimate ( p8| ) above coincides with eq. (|30|), up to a factor 1/V reflecting the difference 
between the real and imaginary parts. 



(Let us parenthetically note that in eq. (22) we substituted the soft cut-off \x\j + k 2 ^ 



by a step function at |fc| = ^o- This is unimportant for the infrared renormalon where 
the integral is saturated at k 2 ~ A*o/ n - However, with the soft cut-off restored the very 
same integral fl22] ) produces the ultraviolet renormalon due to the domain k 2 ~ fj^n. It is 
not difficult to check that the corresponding factorial behavior is the same, up to a sign, 
(C n )uv = (— 2) n n!. In contrast to the situation with the polarization operator IT considered 
in Sect. 2, these two singularities in the Borel plane, infrared and ultraviolet, are symmetric 
with respect to the origin.) 

The statement above - the impossibility of defining nig 16 to the accuracy better than 
Aqcd - implies some tacit assumptions. In particular, we assumed that one should use the 
value of the running coupling a s (k 2 ) in the integrand in eq. ( f22|) . This natural prescription is 
easily justified in QED where the Ward identity reduces the renormalization of the coupling 
constant to the corrections to the photon propagator. In non-abelian theories this is not 
the case in covariant gauges. A general argument below illustrates the fact that one cannot 
get anything else. 

Being interested in effects occurring at the scale /i, much below the mass of the heavy 
quark, one integrates out all momenta above ~ li and arrives at an effective field theory 
of a nonrelativistic heavy quark with an ultraviolet cutoff ll. The parameters of QCD 
are then a s (/i) and ttiq^) (and, in principle, masses for the "light" quarks); including 
external interactions adds also the corresponding couplings that must be renormalized at 
the same scale ll as well. The dependence of these parameters defining the effective theory 
on the renormalization point ll follows from the requirement that physical observables do 
not depend on ll. If the renormalization procedure is such that lowering li from a value ll\ 
down to fi2 incorporates radiative corrections due to virtual momenta \k\ between ll\ and 
Li2, then there must be a linear dependence of mg(/j) on /u, necessarily of the form 

= P m (a s (LL)) = -0$a.(jj.) + ... . (31) 

The fact that it originates from radiative corrections is indicated by the explicit factor a s 
on the right-hand side. Of course, the exact form of the function (3 m (a s ) depends on the 
particular renormalization scheme. If one uses a scheme that coincides to one loop with the 
prescription of introducing a "gluon mass" /U, then according to eq. (TLq) 



= | (32) 



The renormalization group equation (|3l| ) combined with Eq. ( |32|) is equivalent to Eq. (22): 
0(a 2 ) terms in the function j3 m (31) are neglected as subleading. 



pole 



Recalling the discussion in Sect. 2 one is tempted to relate the infrared part in tuq 
in eq. (28) to the matrix element of some local gauge invariant operator. Alas, such an 
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attempt is doomed to fail. This is most easily seen by inspecting the local operators of the 
relevant dimension. The only potential candidate is 



QiDoQ 

where iDq = ido + g s Ao ( we imply gauge invariance plus a static description of the field Q 
similar to that used in HQET). However, the equations of motion reduce this operator to 
those of dimension 5 which can generate corrections of the relative weight h^ CD /rnQ only, 
rather than Aqcd /mQ, provided that the definition of the quark mass is properly adjusted, 



so that there is no so-called residual mass [23], see below. 



From the derivation given above (see Sect. 3 and, especially, the fact that ko = 0, as 
emphasized there) it is clear why the standard OPE program is inapplicable to fn^ le . By 
analyzing Fig. 3 we have realized that only very small frequencies (of the order of {Mq/uiq) 
contribute to the pole mass. In other words, even though the characteristic spatial distances 
are small in the problem at hand, the time separation is parametrically large! To state it 
in more physical terms: measuring the pole mass of the heavy quark requires a very long 
time, inversely proportional to the allowed uncertainty in the absolute value of the mass. 

This means that the non-perturbative infrared contribution in m^ le cannot be expressed 
in terms of a local condensate, but, rather, through a non-local expectation value. We 



have encountered with a similar situation previously [24| in connection with the so-called 



temporal distribution function defined through the hadronic matrix element of the operator 

(DiQ(t,x = 0))e" i /o Ao ^ dT (D t Q(t = 0,3=0)) (33) 

where Q is the heavy quark field and Di is the spatial component of the covariant derivative. 
The pole mass formally appears in consideration of the operator ( |33|) if one considers the 
matrix element of this operator over the heavy quark state in the limit t — > oo. Thus, we 
deal here with a generalization of the standard Wilson path operator for an open path along 
the t direction. 



5 Irrelevance of the Pole Mass 

The only systematic approach presently known that allows one to treat non-perturbative 
effects in QCD analytically is based on Wilson's OPE. This procedure requires - as discussed 
in detail in Sect. 2 - a careful separation of contributions from large and short distance 
dynamics. Our findings from the previous section suggest that the pole mass cannot be 
handled within such an approach. In this section we will demonstrate that this is indeed 
the case. 

To be specific let us discuss the problem of charmless semileptonic decays of B mesons 
which has been already mentioned in Sect. 1. In the parton model the total width is given 
by the probability of the b — > ulPi transition; neglecting all corrections we have 

r. = ^lv-P. (34) 

When a s corrections are included the result of the explicit perturbative computation is 
most naturally expressed in terms of the pole mass (and this is what is usually done, see 
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e.g. ref. @), 



r peri = 1 - -—) + ...). (35) 



This formula is perfectly legitimate as long as we stay in perturbation theory. If, however, 
we would like to take into account non-perturbative (infrared) effects the use of eq. (|35|) 
can be seriously misleading. Let us clarify this point. 

As explained above, the perturbative series for m^ ole diverges factorially - a fact sig- 
nalling the presence of a non-perturbative contribution. A similar factorial divergence takes 
place in the a s expansion in the brackets of eq. (|35|). This perturbative series does have 
the renormalon singularity at 6 = 2ir/b giving rise to (uncontrollable) corrections of or- 
der Aqcd/i^q- Both effects combine, however, to cancel each other, if the running mass 
mQ(mg) is used. 

To see how this cancellation works it is sufficient to consider one-gluon exchange graphs, 
Fig. 5. One must keep in mind that the gluon Green function is assumed to be dressed as 
in Fig. 4, so that a s (k 2 )/k 2 must be used for the gluon propagator. This gives rise to the 
usual infrared renormalon. 

We remind the reader that we are interested in the imaginary parts of the diagrams 
in Fig. 5 corresponding to the appropriate cuts (i.e. the total semileptonic width). The 
graphs 5a and b contain the effect of the radiative shift of the mass discussed in Sect. 3, 
and, in particular, the factorial divergence, see Sect. 4. Our point is that the integrands in 
the Feynman integrals for the diagrams 5a and d will completely compensate each other in 
the domain of virtual momenta \k\ -C rrif,; likewise with the graphs 56 and e. The easiest 
way to demonstrate this compensation is to contract the fermion loop in diagram 5a to a 
local effective vertex where p is the momentum carried by the heavy quark line. The net 
effect of the diagrams 5a and 56 is then to convert the bare 6 quark mass into (iri^ ole ) 5 in 
the total width where m^ ole in this approximation is 

m poie _ m (o) _ £, ^g-j 
and — E is given by the graph of Fig. 3. Explicitly these graphs produce 



Ar 



a+b 



The result for E is implicitly given in eq. (|l6|); its explicit expression is not needed for our 
purposes here, where we have to consider only the integrand of the integral determining E. 

Let us now turn to the diagrams 5d and 5e. First we observe that the vertex we had 
as an effective vertex representing the light fermion loop in Fig. 5a is actually a 

b(iffib 

term in the effective Lagrangian. The gauge invariance of the theory implies that the 
emission of the soft (i.e. \k\ <C mj) gluon from this fermion loop - Figs. 5 d and e - is 
completely taken care of by the substitution 

b(iffib^ b{i]J)fb 
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where is the covariant derivative. Diagrammatically these effective vertices are presented 
in Fig. 6. The gluon emission vertex in Fig. 66 is 

5/4^ 5(m[ 0) ) 4 4. 

After these remarks it is straightforward to see that the diagrams 5d and 5e yield 

AF rf+e = ^^5(m<°>) 4 E. (38) 

We conclude that the uncertain contribution from the domain \k\ -C mb present in individual 
graphs and responsible for the factorial growth of the coefficients (see Sect. 4) is absent in 



the sum, eqs. (37) and (|3q). Of course, the cancellation described above takes place only 
to leading order in k/mb- The residual difference between the graphs 5a and d shows up at 
the level of integrals of the type 

d 3 k/\k\ 



which are harmless from the point of view of the linear in 1/mf, effect we focus on here. Let 
us also note in passing that the diagrams 5c and 5/ do not produce contributions of this 
type and this type of divergence. 

Thus, in applying heavy quark theory one has to avoid the pole mass altogether and to 
use, instead, the running mass which naturally appears in OPE; for only in this case one 
may hope to get consistent and well-defined expansions in powers of 1/mQ. The leading 
operator in the expansion in the problem at hand is b{i Jp^b. From the consideration above 
it follows that the normalization point [i = mj, is the most natural choice: by adopting 
this normalization point we avoid any large logarithms as well as the problem of a factorial 
divergence of the type discussed in Sect. 4. Non-perturbative effects enter through the 
matrix element of this operator; they are also represented by matrix elements of other 
(subleading) operators, for instance, b(i JZ)) 3 iaGb. 

Using the equations of motion one reduces the leading operator to mfbb, where both 
mft and bb are taken at /j = m^. We then evolve bb down to a low normalization point, 
fj. <C mj; the net effect of this evolution is reflected in a factor of the type c(fj,,mb) = 
1 + ai(mb)a s (mb) + a2(mb)a 2 (mb) + .... which is, anyway, included in the perturbative 
calculation. This factor contains also terms of order (///m^)™ due to the exclusion of the 
domain below /j, from the perturbative calculation. It is important that the power n starts 
from n = 2. These small terms, (fj,/mb) n , will contain the series in a s (fi), not a s (mb). 
The additional fi dependence occurring in this way will be canceled by that coming from 
appropriate local operators, e.g. b(i fi)) 3 iaGb. 

Once the operator bb is evolved down to a low normalization point we can use the 
relation 

bb = v^b + -^ biaGb + b((iD) 2 - {ivD) 2 )b + 0(\jm\) (39) 

explicitly demonstrating the absence of 1/mQ corrections. Eq. (^) is valid up to terms 
representing total derivatives. 

In other words, one integrates out the momenta above the scale /j to get a generic 
operator product expansion for the width 

f = c(/i, m b ) ■ mfbblfj, + c 2 (/j, m b ) ■ mqbiaGbl^ + . . . ; (40) 
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the matrix element is taken then by using eq. (|39|), 

A question which immediately comes to one's mind is as follows. What happens if we 
first evolve the operator b(i ]ft) 5 b down to a low normalization point and only then use the 
equations of motion. At first sight we would get mjj (/■*), not m \ (jnb) m this case, and in the 
limit [i — > we would recover (m^°' e ) 5 . The loophole in this argument is rather obvious: 
the operator b{i p) 5 b mixes under renormalization with the operators b(i P) 4 b, etc., and 
accounting for this mixing returns us to the mass normalized at m&. 

It is instructive to dwell on this issue of mixing in more detail, the more so since it is 
intimately related to the notion of the residual mass, to be discussed below. Consider the 
operator hi Pb\^ Q and evolve it down to /u. The effect of the evolution is described by the 
diagrams of Fig. 66. If only terms linear in fj, are kept, the only relevant graph is the vertex 
renormalization; this graph is the same as that for the mass renormalization (Fig. 3), up 
to a sign. As a result we get 



2tt ol 

Upb\n = tiipb\p-— — (w-»jbb. (41) 

O 7T 



Using now the equations of motion in combination with eq. (17) we see that both the left 



and the right-hand sides contain m^/io), q.e.d. A similar relation holds, of course, for any 
power of i Ifi. 

Above we have demonstrated our assertions considering explicitly a certain class of 
diagrams. The result is more general, of course. It can be viewed as a statement that 
the inclusive widths of heavy flavors do belong to the class of observables which are given 
by operator product expansions in the standard understanding: physics of short and large 
distances f^] can be separated into operators and their coefficients, and the infrared behavior 
of the latter is governed, in turn, by the corresponding local operators. The existence of 
the 1/rriQ renormalon in the heavy quark mass shows that this statement is not as trivial 
as it might seem at first sight. However, as soon as the validity of the OPE is accepted, one 
necessarily arrives at the irrelevance of the pole mass. 

Within the OPE-based approach one obtains the corresponding inclusive width as a 
series of operators of the form b...b initially normalized at the scale m^; the infrared stability 
of the inclusive widths ensures that the coefficient functions are finite in any perturbative 
order. Some of these operators contain covariant derivatives acting on the b fields; due to 
the equations of motions these derivatives reduce to powers of the high scale mass mb(rrib). 
Subleading operators can produce relative effects not larger than 1/m^. 

Perturbative corrections actually drop out altogether from certain observable quantities. 
Differences in the lifetimes of different species of hadrons in the same heavy flavor family 
provide a prominent example: the widths of pseudoscalar mesons Pq and baryons Aq agree 
through order 1/mQ 



r , /° ~ l/m 2 Q , 



in clear contrast to their masses, 



Ma q - M Pq 
Ma + M Po 



1/mg. 



3 It must be noted that it is indeed a statement of separation of large and short distances in the process, 
not of perturbative versus, non-perturbative effects; for the Wilson coefficients in general include non- 
perturbative contributions generated for example by instantons of small size ~ 1/fi or 1/rriQ. 
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We plan to present a more detailed discussion in a forthcoming paper J25fl . 

Since our conclusions obviously do not depend on the particular decay process, they 
apply directly and equally to radiative and nonleptonic decays. 



6 The Running of A and other Conclusions 

We have discussed in some detail the problem of the pole mass in the heavy quark the- 
ory. Due to the peculiarities of the static limit an infrared term linear in Aqcd/i^q is 
generated in the pole mass, as signalled by an infrared renormalon. The presence of this 
non-perturbative term makes the notion of the pole mass, beyond perturbation theory, not 
only useless but, rather, detrimental. What is even worse, this non-perturbative term can- 
not be absorbed into any local condensate, unlike the usual OPE-based prescription where 
the infrared effects are naturally incorporated through the condensates. Thus, the pole 
mass should be avoided altogether in analyzing calculable observable quantities. The prob- 
lem disappears provided one uses the running mass thq (//) normalized at a sufficiently high 
point. The same situation in a different context has been noted recently in Ref. [26|. 



The occurrence of a new, linear infrared effect in the static limit reminds us of high- 
temperature QCD. In this theory an external parameter, temperature, sets the scale of the 
energy transfers, and if T ^ ^-qcd this scale is fixed to be of order T. The integrals over 
the four-dimensional momenta degenerate into integrals over three-dimensional momenta. 
This produces new linear infrared divergences in high-temperature QCD [p7f| . Whether this 
parallel leads to non-trivial insights into our problem remains to be seen. 

The fact that corrections linear in 1/ rriQ are present in some quantities is not surprising 
by itself. Perhaps, the best-known example is the axial constant /g (/o or for charm 



and beauty, respectively). In the limit tjiq — > oo they scale with tuq as follows [£8[ 29, 



fB 

f D m b a s (m c ) 

This scaling is known to be significantly modified by pre-asymptotic 1/rriQ terms if one 
defines the axial constant in the standard way via matrix elements of the axial current, 

(0\Ql^q\P Q ) = ifQPv (43) 
However the correction is much smaller if one would define it via the pseudoscalar current 



32fl: 



(0\Qi^q\P Q ) = IqM Pq (44) 
In the heavy quark limit the two definitions obviously coincide, 

k ^ ggg (45) 
Jq m Q + m q 

and this difference is therefore indeed contained in 1/m.Q terms. 

Calculation of /q cannot be directly formulated as an OPE-based procedure; therefore, 
the emergence of the 1/mn correction is natural. The pole mass belongs to the same class. 
Yet the OPE prescription is fully applicable to the calculation of the inclusive widths of 
heavy flavor hadrons; the infrared behavior of the corresponding Wilson coefficients is given 
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by matrix elements of the appropriate operators. This difference between the pole mass 
and the inclusive widths is not accidental, of course, and could be anticipated. 

A remark is in order here concerning the place the pole mass occupies in the heavy 
quark effective theory |3j, ||| . HQET is formulated in such a way as to get rid of the heavy 
quark mass at all; the effective theory is left without a large parameter which might set 
an appropriate scale for distances. It is then only the renormalization point /x that fixes 
the scale of momenta for static quantities. However the pole mass assumes that one takes 
the limit fj, — > and, thus, no parameter is left at all. It is clear that in any consistent 
formulation of HQET it is impossible to set /i = 0. One should keep /i explicitly and operate 
only with mg( / u). 

It is well known that in the framework of the HQET a key role is played by the difference 
between the the hadron and the quark mass for an asymptotically heavy quark: 

X = 4^J^--Q) • (46) 

It is always stated that the mass tjiq entering this definition is the pole mass. We have 
shown, however, that this quantity is ill-defined. There have been a few attempts to define 
it consistently in HQET. One of those has been made in ref. | |23| |; namely, for pseudoscalar 
mesons it was suggested that 

_ {0\i(vd)(gi l5 h v )\P(v)) 
P (0\qij 5 h v \P(v)) ' 1 > 

The fields of the HQET are assumed here and standard notations are used. This def- 
inition per se is not better and is plagued by the same problems - the impossibility of 
disentangling non-perturbative effects from the perturbative contributions. Any consistent 
formulation has to discriminate large and small momenta rather than perturbative versus 
non-perturbative effects. A consistent formulation must then include the renormalization 
point /x to ensure that momenta higher than [i do not appear in the matrix elements. If /i 
is introduced one must then define a "running" value of A(/i), depending on the renormal- 
ization point /z, as follows 

A(/i) = J™JM Hq - m Q (/i)) . (48) 

with fi having to exceed sufficiently typical hadronic scales. Its renormalization point de- 
pendence is given by eq. (31) ; attempts to put /i to zero to arrive at the 'old' A, would 
bring about uncontrollable uncertainties of order Aqcd in A(0). 

It has been suggested [ 33 1 to use the requirement that the so-called residual mass term 
vanishes to rigorously define the heavy quark pole mass rn^ le . In Sect. 5 it was shown that 
mixing between the operators Qi I/)Q and QQ arises already at the one-loop level to order 
[i. This means that even if the effective Lagrangian of HQET is chosen in such a way that 
at a certain /i the residual mass is zero, it necessarily re-appears at a different value of /i. 
If the process is characterized by a single scale (like txiq in the total inclusive widths) one 
can certainly adjust the effective Lagrangian so that there is no residual mass. If the scale 
varies, however, it is mandatory to "readjust" the notion of the mechanical mass of the 
heavy quark. In other words by requiring the vanishing of the residual mass one defines the 
running quark mass. (It has been traced [34] how various observables calculated in HQET 
via matching with full QCD turn out to be independent of this term although it is present 
in intermediate stages Q) 

4 N.U. is grateful to V.M. Braun for pointing out this constructive interpretation. 
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The fact that perturbative corrections to mg vanish in dimensional regularization, 
which is often referred to as a remedy, is an accident due to the fact that HQET does 
not contain a dimensional parameter at the perturbative level. References to dimen- 
sional regularization per se without an actual subtraction scheme is irrelevant for the ef- 
fective theory which requires a clear separation of high and low momenta. For example, 
the MS scheme giving a logarithmically dependent one-loop value for niQ in the form 
niQ(n) = mQ(mg)(l + (2a s /7r) log(mg/ju)) for an infrared cut-off fi <C mg cannot achieve 
the momentum separation necessary for OPE. 

Note that the problem of the total inclusive widths is, strictly speaking, outside pure 
HQET: for the large mass scale parameter - the energy release - is intrinsically involved 
from the very beginning, determining the characteristic space-time separations here. 

Some relations are known for exclusive formfactors in b — > c transitions, that, in the 
limit mb,m c — > oo depend on A outside the zero recoil kinematics (see e.g. p5|| ). However 
they contain purely perturbative 0(ct™) corrections. Their main piece reflects the hybrid 
renormalization |3(], [3(| that must be properly accounted for in constructing the effective 
low energy description; they are governed by the scale momenta between nib an d m c , or, in 
general, above /i. Some corrections, appearing already in the effective theory itself, reflect, 
however, much lower scales; obviously in the heavy quark limit at v ^ v' there always 
exists a large domain of gluon momenta ^ \k\ <J |w|ttiq whose effect is not governed by 
a s (mg), but rather by the running coupling at the lower scale. These corrections therefore 
can lead again to infrared renormalon contributions. There is no reason for them not to 
be able to convert the "pole" A into a 'running' one. It is clear that this problem deserves 
further theoretical studies, especially if one wants to understand the real meaning of the 
relations mentioned above in the presence of an infrared renormalon in the heavy quark 
mass. In principle a more complicated situation is conceivable - exclusive modes may not 
be described completely by the standard OPE and intrinsic uncertainties similar to that in 
the pole mass could have emerged. 

We do not address in this paper the numerical aspects of the infrared renormalons; they 
are left for the future studies |^5| (see also Still for orientation it is instructive to 



consider the estimate (28). According to this expression the uncertainty in the pole mass is 



5m p Q ole ~ ^A QCD . (49) 

To determine what particular A (i.e. A corresponding to what particular subtraction 
scheme) enters on the right-hand side one needs to perform a two-loop calculation. In 
the absence of such calculations it is consistent to use the one-loop A which is one and the 
same in all schemes, A one i oop ~ 100 to 150 MeV. Then one gets 

5m p ° le ~ 50MeV; (50) 

certainly this estimate needs improvement via real two loop calculations. 

The presence of the l/rriQ renormalon in the pole mass of the heavy quark makes the 
attempts to extract the value of the pole mass from the absolute widths of heavy flavors 
[ j38| , [39] , |Iof not very meaningful from a theoretical perspective; it is the masses m c {w%) and 
mb(ml) (see ref. p8| ) that can be extracted from the widths. Indeed, one can reach the 
required 1/mQ accuracy through calculating sufficiently many perturbative corrections only 
if the semileptonic width of D or B mesons is expressed in terms of the high scale mass. On 
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the other hand, one could determine V c b from the b — ► c semileptonic width assuming that 
one is close to the small velocity (SV) limit @ - then this width depends mainly on the 
difference m\, — m c meaning that the overall uncertainty in the mass cancels out in m& — m c 



The uncertainty in the pole mass can well constitute a sizable part of the commonly 
accepted value for A for mesons of about 400 -f- 600 MeV, cf. eq. (|50"D . This is apparently 
due to the peculiar Coulomb enhancement of the leading radiative corrections that has been 
noted in a different context by Braun et al. [41]. 
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Figure Captions 

Fig. 1 . Feynman diagrams for U^ u in eq. (2). Solid lines denote the quark q, dashed 
lines gluons. 

Fig. 2 . "Bubbles" in the gluon Green function giving rise to the infrared renormalon 

in 

Fig. 3 . One-loop diagram for the mass renormalization. Thick line denotes the heavy 
quark. 

Fig. 4 . Infrared renormalon in the pole mass. 

Fig. 5 . Diagrams determining 0(a s ) corrections to the total semileptonic width (taking 
the imaginary part is implied). 

Fig. 6 . Graphs 5a and d in the language of the effective local vertices (denoted by the 
black box). 
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